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Abstract 

We study three generation models in the four-dimensional spacetime, which can be derived 
from the ten-dimensional N = 1 super Yang-Mills theory on the orbifold background with 
a non-vanishing magnetic flux. We classify the flavor structures and show possible patterns 
of Yukawa matrices. Some examples of numerical studies are also shown. 



* E-mail address 
^E-mail address 
"'■E-mail address 
^E-mail address 



abe@ tuhcp . phys . tohoku .ac.jp 
kschoi@gauge . scphys . kyoto-u .ac.jp 
kobayash@gauge . scphys . kyoto-u. ac . jp 
ohki@scphys .kyoto-u .ac.jp 



1 Introduction 



Extra dimensional field theories, in particular string-derived extra dimensional field theories, 
play important roles in particle physics, e.g. as an origin of the flavor structure including 
the hierarchy of quark/lepton masses and mixing angles. How to derive chiral theory is a 
key issue when our starting point is extra dimensional theory. Introduction of a magnetic 
flux in extra dimensional space is one of interesting ways to obtain chiral theory. Indeed, 
several studies have been carried out on models with magnetic fluxes in field theories and 
superstring theories [HEIEIIIIEIEIITIIH]. Furthermore, magnetized D-brane models are T- 
duals of intersecting D-brane models and within the latter framework several interesting models 
have been constructed [2 El El El [TOl fTlfl . 

Zero-modes are quasi-localized on the torus with the magnetic flux. The number of zero- 
modes, which corresponds to the generation number, is determined by the value of the magnetic 
flux in the same way as that the generation number is determined by the intersecting number in 
intersecting D-brane models. Yukawa couplings among zero-modes in four-dimensional effective 
theory are obtained by overlap integral of zero-mode profiles in extra dimensions. Large or 
suppressed Yukawa couplings can be derived depending on the size of overlap integral. That is, 
when zero-modes are quasi-localized far away from each other, their couplings in 4D effective 
field theory are suppressed. On the other hand, when their localized points are close to each 
other, 4D effective Yukawa couplings would be of 0{1). Thus, magnetized torus models are 
quite interesting to derive realistic models, in particular a realistic flavor structure. However, 
it is still a challenging issue to derive realistic mass matrices of quarks and leptons. 

Orbifolding the extra dimensions is another way to derive chiral theory pTSj. In Ref. [H], 
magnetized orbifold models have been studiecl^. Phenomenological aspects in magnetized orb- 
ifold models are different from those in magnetized torus models. Some of zero-modes are 
projected out by the orbifold projection. However, odd modes as well as even modes could 
correspond to zero-modes, although odd modes correspond to only massive modes on orbifolds 
without the magnetic flux. Then, the generation number is smaller than the number of the 
magnetic flux, i.e. one in magnetized torus models with the same magnetic flux. Thus, a new 
type of flavor structure can appear in magnetized orbifold models. Hence, it is quite impor- 
tant to study in detail phenomenological aspects of magnetized orbifold models. That is our 
purpose in this paper. We classify three generation models and study predicted patterns of 
Yukawa matrices. 

The paper is organized as follows. In section 2, we give a review on magnetized orbifold 
models. In section 3, we classify three generation models on the orbifold with magnetic fluxes. 
In section 4, we study Yukawa couplings in three generation models and we show explicitly an 
example of numerical studies on our models. Section 5 is devoted to conclusion and discussion. 
In Appendix, we show explicitly all of possible Yukawa matrices in our three generation models. 

^ See for a review [12] and references therein. 

^ Other geometrical backgrounds with a magnetic flux have also been studied [TSl I16j. 



2 Magnetized extra dimensions 



Here, we give a review on extra dimensional models with a magnetic flux on torus and orbifold 
backgrounds [3 [H] . 

2.1 U{N) gauge theory on (r^)^ 

We start with A/" = 1 ten-dimensional U (N) super Yang-Mills theory. We consider the back- 
ground R^'^ X (T^)'^, whose coordinates are denoted by (/i = 0, ■ ■ ■ , 3) for the uncompact 
space R^'^ and i/m {rn = 4, ■ ■ ■ ,9) for the compact space (T^)^. At the first stage, we use orthog- 
onal coordinates of the compact space and choose the torus metric such that i/m is identified by 
Um + nm with rim = integer, i.e. Um ^ Vm + ^- At the end of this subsection, we will extend it 
by introducing the complex structure. Also, we can extend the following discussions to A/" = 1 
super Yang-Mills theory on i?^'^ x (T^)". 
The Lagrangian is given by 



where M, = 0, ■ ■ ■ , 9. Here, A denotes gaugino fields, F*^ is the gamma matrix for ten- 
dimensions and the covariant derivative Dm is given as 

Dm\ = OmX-iIAmAI (1) 

where Am is the vector field. Furthermore, the field strength Fmn is given by 

Fmn = OmAn - On Am - i[AM, An]. (2) 

The gaugino fields A and the vector fields A^ corresponding to the compact directions are 
decomposed as 

Kx,y) = ^Xn{,x) ®il)n{,y), 

n 

Am{x,y) = '^iPn,m{x) ® (f)n,miy)- 
n 

Here, we concentrate on zero-modes, '?/'o(z/) and we denote them as '?/'(?/) by omitting the sub- 
script "0" . Furthermore, the internal part ipi^y) is decomposed as a product of the i-th parts, 
i.e. ip(i){y2i+2,y2i+3)- Each of V^(i)(?/2i+2, 1/21+3) is two-component spinor, = (^(j)4^, 
and their chirality for the i-th part is denoted by Sj. We use the gamma matrix F'" corre- 
sponding to the i-th as 

^ -[10)' ^ - [t ^"^^ 

and the total gamma matrices are obtained as their direct products with the four-dimensional 
part. 
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Here, we introduce the magnetic flux in the background as F45, Fgy and Fgg, which are given 



by 







^45 



^67 



89 



27r 



27r 



27r 



ATixAfi 







V 







V 







This background breaks the gauge group U{N) as U{N) YZ=iU{Na) with = Ea^a- We 
concentrate on an Abehan flux, although in general non-Abelian magnetic fluxes, which reduce 
ranks of gauge groups, are possible |T7l [HI [19] . 

Here we focus on the U{Na) x U{Nb) part, which has the magnetic flux. 



-p2i+2,2i+3 — 27r 











for i = 1, 2, 3. We use the following gauge, 

^2i+2 = 0, ^2i+3 = -^2j+2,2i+3 Z/2i+2- 

Similarly, the gaugino fields A and their i-th torus parts are decomposed as 



X{x,y) 



(4) 



(5) 



(6) 



The fields A"" and A correspond to the gaugino fields under the unbroken gauge group U{Na) x 
U{Nb). On the other hand. A"* and A'"^ correspond to bi-fundamental matter fields, {Na,Nb) 
and {Na,Nb), under the unbroken gauge group U{Na) x U{Nb). The Dirac equations for these 
gaugino fields corresponding to zero-modes are obtained as 



diip 



aa 
(04 



[a, + 2vr(MW-M«)l/2.+2]V^[^V ^ 



[a. + 27r(M«-M»)y2.+2]V^: 



ba 

(i)H 
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(i)H 



(i)- 



[a,-27r(M«-M«)y2.+2]^[^V ^ 



1^ [9,-27r(M«-M«)y2.+2]^ 



ba 

(i)- 



diip 



bb 
(i)- 
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where Bi = 821+2 + id2i+z and di = 821+2 — id2i+3- The gaugino fields, ip"'"' and ip^^, for the 
unbroken gauge symmetry have no effect from the magnetic flux in their Dirac equations. 
Hence, they have the same zero-modes as those on (T^)^ without the magnetic flux. On the 
other hand, the magnetic flux appears in the zero-mode equations of ip"'^ and ip^"" corresponding 
to bi-fundamental matter flelds, {Na, Nb) and {Na, Nb). Furthermore, they satisfy the following 
boundary conditions, 

= Cf(2/2i+2,2/2i+3), 
= ^J(l/2i+2,l/2i+3), 

because of Eq. ([5]). 

For the i-th with M^^ -M^^ > 0, the fields ^/^g+ and have |M^*) -M^'^ \ normahzable 
zero-modes, while ^'(f)- and 4'(i)+ have no normahzable zero-modes. Thus, we can derive chiral 

theory. When M^*^ — Mjf '' < 0, ipf^y and 'ip'^i^_^_ have |M^*^ — M^^*^| normahzable zero-modes. 
The normahzable wavefunction for the j-th zero mode is obtained as 



C(?/2i+2 + l,y2i+3) 

^?fe+2 + l,y2i+3) 
^s-{y2i+2,y2i+3 + l] 
^sf(l/2i+2,Z/2i+3 + 1) 



0^(l/2i+2,l/2i+3) 



-Mtv 



y2i+2^ 



j/M 




iM{y2i+2+iy2i+3),Ml) 



for M = |M« - M^^ I and j = 0, 1 



— 1, where Nj is a normalization constant and 



j/M 




{M{y2i+2+ty2i+3),Mi) 



-MTT{n+j /Mf+2wi(n+j /M)M{y2i+2+iy2i+3) 



that is, the Jacobi theta-function. Furthermore, we can introduce the complex structure mod- 
ulus r by replacing the above Jacobi theta-function as 



j/M 




{M{y2^+2 + iy2i+3),Mi) 



j/M 




iM{y2i+2 + Ty2^+3),MT). 



The total number of bi-fundamental zero- modes is given by nLi l^i'^ -^b'^ I and all of them 



Since the ten-dimensional 



have the same six-dimensional chirality sign nf=i (^i*^ ~ 
chirality of gaugino fields is fixed, bi-fundamental zero-modes for either {Na,Nb) ot {Na,Nb) 
appear for a fixed four- dimensional chirality. That is, the total number of bi-fundamental 
zero-modes for {Na, Nb) is equal to 

3 
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i=l 



When lab < 0, this means that there appear \Iab\ independent zero modes for {Na,Nb). It is 

M(i) - Mj:'\ Zero -mode wavefunctions are 



also convenient to introduce the notation, J*^ 
given by a product of two-dimensional parts, i.e. 
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for zi = 0,---,(|M« -M^'^l -1), Z2 = 0,---,(|M(2) -Mf I - 1) and zs = 0,---,i\MP - 
Mf)| - 1). 

2.2 U{N) gauge theory on magnetized orbifolds T^/{Z2 x Z2) 
Here we review on the T^/{Z2 x orbifold with a magnetic flux [H] . 

2.2.1 TVZ2 orbifold 

First, let us study the f/(A^) gauge theory on the orbifold with the coordinates (1/4,1/5), 

which transform as 

Vi -z/4, 1/5 ^ -1/5, 

under the Z2 orbifold twist. Here, we associate the Z2 twist with the Z2 action in the gauge 
space as 

A^{x,-y) = PA^{x,y)p-\ Am{x,y) = -PA^{x,y)p-\ 

and the Z2 boundary conditions for gaugino fields, 

A±(x,-i/) = ±PXi{x,y)p-\ 

where the Z2 projection P must satisfy P^ = 1. 

We focus on the U{Na) x U{Ni,) block (g]), ([6]) and consider the spinor fields, A^'^, Xf, A^" 
and A^, in particular bi-fundamental fields A^^ and A^, where ± denotes the chirality Sj in 
the extra dimension. Without the Z2 projection, there are \Ma — Mh\ zero modes for A^* and 
^6a_ Pqj, example, when Ma — Mf, > 0, A"^ as well as A^" has {Ma — Mb) zero modes with the 
wavefunctions B-^ for j = 0, ■ ■ ■ , {Ma — Mb — 1). When we consider the Z2 projection, either 
even or odd modes of them remain. Here note that 

Q\-y,,-y,) = 6*^-^(1/4,1/5), 
where 6*^(1/4,1/5) = 0°(l/4,l/5)- That is, even and odd functions are given by 

©Ln = ^(e^ + e*^-^), 

e^odd = ^(e^-e*^-^), (7) 

respectively. For example, when we consider the projection P such that X'^{x, —y) = X"^{x^y), 
only zero-modes corresponding to 6^ven remain and the number of zero-modes is equal to 
{Ma - Mb)/2 + 1 for {Ma - Mb) = even and {Ma - Mb + l)/2 for {Ma - Mb) = odd. On the 
other hand, when we consider the projection P such that X'^{x, —y) = —X°^{x,y), only zero- 
modes corresponding to Oodd remain and the number of zero- modes is equal to {Ma — Mb)/2 — 1 
for {Ma - Mb) = even and {Ma - Mb - l)/2 for {Ma - Mb) = odd. The same holds true for 
A'L"- Table 1 shows the numbers of zero-modes with even and odd wavefunctions for M < 10. 
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M 


0123456789 10 


even 


1122334455 6 


odd 


0001122334 4 



Table 1: The numbers of zero-modes for even and odd wavef unctions. 
2.2.2 TV(Z2 X Z'^) 

Here, we can extend the previous analysis on the two-dimensional orbifold T"^ /Z2 to the U{N) 
gauge theory on the six- dimensional orbifold T^/ [Z2 x Z'2). We consider two independent twists, 
Z2 and Z'2. The Z2 twist acts on the six-dimensional coordinates Hm (m = 4, ■ ■ ■ , 9) as 

Vm -Vm (for m = 4,5,6, 7), y„ yn (for n = 8,9), 

and the Z2 twist acts as 

Vm -Vm (for m = 4,5,8,9), y„ yn (for n = 6, 7). 

If the magnetic flux is vanishing, we realize four-dimensional Af — 1 supersymmetric gauge 
theories for the orbifold /{Z2 x Z2). The bi-fundamental matter fields ^^1*^2, S3) -^51,52,53 '^i^^i 
the chirality Si corresponding to the i-th are also introduced. Their Z2 boundary conditions 
are given by 

'^si,S2,S3 (-^j ymiyn) — -^51,82, S3 (-^j Z/m; Z/n)-f i 

With m = 4,5,6,7 and n = 8, 9 for A»^%^,,3, A^«,^,,3 and A^^,,^,,3. Similarly, the Z'2 

boundary conditions are given by 

-^si,S2,S3(^) yrmlln) — ^l^sP ^si,S2,S3{^i Vmi VnjP i 

with m — 4, 5, 8, 9 and n — 6, 7. Then, depending on the projections P and P', even 
or odd modes for the i-th torus remain such as ©^^^ or O^'j^. Their products such as 
ni=3©eve^odd(y2i+2,y2i+3) provide with zero-modes on the 7^/(^2 x Z'2). 



3 Three generation magnetized orbifold models 

In this section, we consider the U{Na) x U{Nb) x U{Nc) models, which lead to three families of 
bi-fundamental matter fields, (iVa, N^) and {Ng,, Nc). Such a gauge group is derived by starting 
with the U (N) group and introducing the following form of the magnetic fiux, 

\ MWl^^xiv. J 
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\ 
\ 

where N = Na + Nf, + N^. For A^^^ = 4, A^^^ = 2 and Nc = 2, we can realize the Pati-Salam gauge 
group up to U{1) factors, some of which may be anomalous and become massive by the Green- 
Schwarz mechanism. Then, the bi-fundamental matter fields, (A^^^, Nt) and {Na, -^c) correspond 
to left-handed and right-handed matter fields. In addition, the bi-fundamental matter fields 
{Nb, Nc) correspond to higgsino fields. We assume that supersymmetry is preserved at least 
locally at the a — b sector, b — c sector and c — a sectorjf] Then, the number of Higgs scalar 
fields are the same as the number of higgsino fields. There are no tachyonic modes at the 
tree level. Indeed, in intersecting D-brane models it would be one of convenient ways towards 
realistic models to derive the Pati-Salam model at some stage and to break the gauge group 
to the group SU{3) x SU{2)l x f/(l). (See e.g. Ref. [m [20] and references therein.fl At 
the end of this section, we give a comment on breaking of SU{4:) x SU{2)l x SU{2)r to 
SUi?,) X SU{2)l X f/(l). 

In both cases with and without orbifolding, the total number of chiral matter fields is a 
product of the numbers of zero-modes corresponding to the i-th for i = 1,2, 3. That is, the 
three generations are realized in the models, where the i-th has three zero-modes while each 
of the other tori has a single zero mode. Thus, there are two types of flavor structures. That 
is, in one type the three zero- modes corresponding to both left-handed matter fields {Na, Nb) 
and right-handed matter fields [Na, N^) appear in the same i-th T^, while each of the other tori 
has a single zero-mode for [Na, Nb) as well as {Na, Nc). In the other type, three zero-modes of 
{Na, Nb) and {Na, Nc) are originated from different tori. The Yukawa coupling for 4D effective 
field theory is evaluated by the following overlap integral of zero-mode wavefunctions [22] 

Yij = J d^y'^Li{y)'^Rj{y)<pH{y), 

where 4'L{y), i'R{y) and 4>H{y) denote zero- mode wave-functions of the left-handed, right- 
handed matter fields and Higgs field, respectively. Note that the integral corresponding to each 
torus is factorized in the Yukawa coupling. In the second type of flavor structure, one obtains 
the following form of Yukawa matrices, 

at the tree-level, because the flavor structure of left-handed and right-handed matter fields are 
originated from different tori. This matrix, Yij, has rank one and that is not phenomenologically 

^See for the supersymmetric conditions e.g. Ref. [71[S]. 

^ See for the Pati-Salam model in heterotic orbifold models e.g. Ref. [H], where SU{4) x SU{2)l x SU{2)r 
is broken to the standard gauge group by vacuum expectation values of scalar fields, (4, 1, 2) and (4, 1, 2), while 
in the intersecting D-brane models SU{A) x SU{2)l x SU{2)r is broken by splitting D-branes, that is, vacuum 
expectation values of adjoint scalar fields. 



"67 



89 



27r 



2n 



V 





Mb'^^lNtxNt 



7 









^bc 


I 


even 


even 


even 


II 


even 


odd 


odd 


ir 


odd 


even 


odd 


III 


odd 


odd 


even 



Table 2: Possible patterns of wavef unctions with non- vanishing Yukawa couplings for the first 
torus. 

interesting, unless certain corrections appear. Hence, we concentrate on the first type of the 
flavor structure. In the first type, the flavor structure is originated from the single torus, where 
both three zero-modes of {Na, Nb) and {Na, Nc) appear. We assign this torus with the first torus. 
On the other hand, the other tori, the second and third tori, do not lead to flavor- dependent 
aspects. That is, Yukawa matrices are obtained as the following form, 

where the structure of ajj"* is determined by only the first torus corresponding to three zero- 
modes (Na, Nh) and {Na, N^) while the other tori contribute to overall factors a*-^-* and a^^\ 
Thus, we concentrate on the single torus, where both of three zero- modes {Na, Nb) and {Na, N^) 
appear, i.e. the first torus. 

Zero-mode wavefunctions are classified into even and odd modes under the Z2 twist. Only 
even or odd modes remain through the orbifold projection. Furthermore, the 4D Yukawa 
couplings are non-vanishing for combinations among (even, even, even) wavefunctions and 
(even, odd, odd) wavefunctions, while Yukawa couplings vanish for combinations among (even, 
even, odd) wavefunctions and (odd, odd, odd) wavefunctions. Thus, we study only the former 
case with non-vanishing Yukawa couplings, that is, the combinations among (even, even, even) 
wavefunctions and (even, odd, odd) wavefunctions. Hence, we are interested in four types 
of combinations of wavefunctions for the first torus, as shown in Table [2l The IF type of 
combinations is obtained by exchanging the left and right-handed matter fields in the II type. 
Thus, we study explicitly the three types, I, II and HI. 

We can realize three even zero- modes when = 4,5, as shown in Table [H On the 

other hand, three odd zero-modes can appear when = 7,8. Furthermore, the consistency 
condition on magnetic fluxes requires 

|j(l)| _ |r(l)| i |r(l)| 

Thus, the number of Higgs and higgsino fields are constrained. Table [3] shows all of possible 
magnetic fluxes for the three types, I, II and HI. The fourth and fifth columns of the table show 
possible sizes of magnetic fluxes for {1^1^ \ and the number of zero-modes corresponding to the 
Higgs fields. As a result, flavor structures of our models with Yukawa couplings are classified 
into 20 classes. However, the model with (|/afc^|, \Icci\, \Ibc\) = (5,7,2) has no zero-modes for 
the Higgs fields. Thus, we do not consider this case, but we will study the other 19 classes 
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U 
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u 
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n 
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II 
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11 


5 
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1 9 


U 
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12 


5 




5 
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13 


6 
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1 


III 
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7 


14 


8 




7 


7 





1 




7 


8 


15 


8 




7 

8 
8 


8 
8 
8 


1 

16 



1 
9 
1 



Table 3: The number of Higgs fields of (T^)^ with non-vanishing Yukawa couplings. 

in Table [31 Therefore, we study possible flavor structures explicitly by deriving the coupling 
selection rule and evaluating values of Yukawa couplings in these 19 classes. That is the purpose 
of the next section. 

Before explicit study on flavor structures of 19 classes in the next section, we give a comment 
on breaking of SU{A) x SU{2)l x SU{2)ji. At any rate, we need the SU{3) x SU{2)l x U{1) 
gauge group at low energy. When the magnetic flux and orbifold projections lead to the 
SU{4:) X SU{2)l X SU{2)f( gauge group from f/(8) as we have discussed so far, we need further 
breaking of 5t/(4) x SU{2)l x SU{2)r to SU{3) x SU{2)l x U{1). Such breaking can be realized 
by assuming non-vanishing vacuum expectation values (VEVs) of Higgs fields like adjoint scalar 
fields for SU{4:) and SU{2)ji and/or bi-fundamental scalar fields like (4,1,2) and (4,1,2) on 
fixed points. Note that our models have degree of freedom to add any modes at the fixed points 
from the viewpoint of point particle field theory. The above breaking may affect the structure 
of Yukawa matrices as higher dimensional operators. However, we will show results on Yukawa 
matrices without such corrections. 

Alternatively, magnetic fluxes and/or orbifold projections break U{8) into U{3) x f/(l)i x 
f/(2)L X f/(l)2 X U{1)3. The gauge group f/(3) x f/(l)i would correspond to f/(4) and U{1)2 x 
f/(l)3 would correspond to [7(2)^. We assume that all the bi-fundamental matter fields under 
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f/(3) X i.e. extra colored modes, are projected out. The bi-fundamental matter fields 

for f/(3) X 17(1)2 and U{3) x U^l)^ correspond to up and down sectors of right-handed quarks, 
respectively. Similarly, up and down sectors of Higgs fields and right-handed charged leptons 
and neutrinos are obtained. In this case, the classification of this section and patterns of Yukawa 
matrices, which will be studied in the next section and Appendix, are available for up-sector 
and down-sector quarks as well as the lepton sector. However, the up sector and down sector 
can correspond to different classes of Table [31 On the other hand, the up sector and down 
sector correspond to the same class in Table [3], when the SU^i) x SU{2)l x SU{2)ji is broken 
by VEVs of Higgs fields on fixed points as discussed above. 



4 Yukawa couplings in three generation models 
4.1 Yukawa interactions 

Following [71 [23], first we show computation of Yukawa interactions on the torus with the 
magnetic fiux. Omitting the gauge structure and spinor structure, the Yukawa coupling among 
left, right-handed matter fields and Higgs field corresponding to three zero-mode wavefunctions. 



e 



i,Mi 1 



e^'^^^{z) and {e'''^^^{z))*, is written by 



ijk 



where z = X4 + ryr,, Mi = M2 = M3 = J^^-* and c is a fiavor- independent contribution 
due to the other tori. Note that Mi + M2 = M3. Because of the gauge invariance, not the 
wavefunction ©'^'^^^(2), but (6'^'*^^ (2))* appears in the Yukawa coupling [7]. 
By using the formula of the ?9 function. 



r/Ni 




[zi,NiT) X ^ 

r+s+Nim 





S/N2 





xt9 



N2r-Nis+NiN2m 
NiN2{Ni+N2) 





Z2, N2T) 
{zi + Z2,t{Ni+N2)) 

{Z1N2 - Z2N1, tNiN2{Ni + N2)) 



we can decompose 0*'*^^(z)0-''*^2(z) as 
Q''^'{z)Q^^^^'{z) = 



Wavefunctions satisfy the orthogonal condition 

dzdzQ^'^^e^'^^y 



M2i-M-ij+MiM2m 
M1M2M3 





(0,rMiM2M3 
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Then, the integral of three wavefunctions is represented by 



c 

IM3I-1 

m=0 



M2i-M-Lj+MiM2m 
M1M2M3 





(0,rMiM2M3) X 6i+j+Mim,k+M3i, 



where i = integer. Thus, we have the selection rule for allowed Yukawa couplings as 

i+j = k, 

where i,j and k are defined up to mod Mi, M2 and M3, respectively^ In addition, the Yukawa 
coupling Yijk, in particular its flavor-dependent part, is written by the function. When 
g.c.d.{Mi, M3) = 1, a signle function appears in Yijk- When g.c.d.{Mi, M3) = g ^1, g terms 
appear in Yijk as 



Y 



ijk 



9 
n=l 



M2k-M-ij+M2Mslo 
M1M2M3 





(0, TM1M2M3 



where £q is an integer corresponding to a particular solution of M3/0 = MirriQ + i + j — k with 
integer mo- 

Zero-mode wavefunctions on the orbifold with the magnetic flux are obtained as even or 
odd linear combinations of wavefunctions on the torus with the magnetic flux ([7]). Thus, it is 
straightforward to extend the above computations of Yukawa couplings on the torus to Yukawa 
couplings on the orbifold. As a result, Yukawa couplings on the orbifold are obtained as proper 
linear combinations of Yukawa couplings on the torus, i.e. linear combinations of functions. 
Here we introduce the following short notation for the Yukawa coupling. 



where 



Vn 



M 



TV 
M 





(0, tM) 



(9) 



M1M2M3. 



Since the value of M is unique in one model, we omit the value of M as well as r for a compact 
presentation of long equations. 

Four models in Table [3] has = 0, where the Higgs zero-mode corresponds to the even 
function, that is, the constant profile. We can repeat the above calculation for this case, that is, 
the case where, one of wavefunctions in ([8]), e.g. Q'^'^^{z) is constant. As a result, the Yukawa 
matrix is proportional to the (3 x 3) unit matrix, Yjk = c'6jk- That is not realistic. Thus, we 
will not consider such models. 

At any rate, we can apply the above selection rule and r/jv for 20 classes of models, which 
have been classified in section 3, in order to analyze explicitly all of possible patterns of Yukawa 
matrices. In the next subsection, we show one example of Yukawa matrix among 20 classes of 
models. In Appendix, we show all of possible Yukawa matrices for 15 classes of models in Table 
[3] except models with I^]^ = and the model without zero-modes for the Higgs fields. 



See for the selection rule in intersecting D-brane models, e.g. Ref. [24l I25j. 
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i?,(A™) 





i 
2 
3 
4 
5 
6 
7 


^ (B ' - B • ) 


^ (BV^ - B«.0 

^ (B ' - B • ) 
(03.7 _ 04,7) 


00,14 

1 tOilM 1 C\13,14\ 
(0 ' + ^' j 

^ (02.14 + 012,14) 
_!= (03.14 ^ 011,14) 
_1= (04.14 + 010,14) 
^ (05.14 + 09,14) 
_1= (06.14 + 08,14) 
07,14 



Table 4: Zero-mode wavefunctions in the 7-7-14 model. 



4.2 An illustrating example: 7-7-14 model 

Let us study the model with \Ica'\, \Ibc^\) = (7, 7, 14). Following Tabled we consider the 

combination of zero-mode wavefunctions, where zero-modes of left and right-handed matter 
fields and Higgs fields correspond to odd, odd and even wavefunctions, respectively. Their 
wavefunctions are shown in Table |H Hereafter, for concreteness, we denote left and right- 
handed matter fields and Higgs fields by Li, Rj and Hk, respectively. This model has eight 
zero-modes for Higgs fields. 

Then, their Yukawa couplings YijkLiRjHk are written by 



where 



ijk-^-^ k 



4 



( -Vc 

- 

V 





72^9 






- 




72^f 




ylj 



-pi- \ 










' V2^e 



'73^^ 



/ 


73^" 





"73^'^ 




-73^^ 




\ 



73^/ 



73^/ 



/ 



T2y^ 
\ 73^/ 



73^^ 
73^^ 73^/ ^ 

73^^^ 





(10) 



and 



Va = V0 + 2?798 + 2?7i96 + 27/294, 
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Figure 1: The iV-dependence of log;^?77v in the 7-7-14 model (M = 686), where A = 0.22 is 
chosen to the Cabibbo angle. The solid, dashed and dotted curves correspond to r = i, 1.5i 
and 0.5i, respectively. Note that tin has a periodicity rj^+nM = Vn with an integer n. 

Ub = ?77 + ^791 + ^7105 + ^7189 + ^7203 + ^?287 + ^7301, 

Vc = ^714 + ^?84 + ^?112 + ^7182 + ^?210 + ^7280 + ^?308, 

Vd = ^21 + ^777 + ^7119 + ^175 + ^7217 + ^273 + ^7315, 

Ve = ?728 + ^?70 + ^?126 + ^7168 + ^?224 + ?7266 + ^?322, 

Vf = ^735 +^63 +^?133 + ^7161 +^7259 + to, 

Vg = ^?42 + V56 + VlAO + ^154 + ?7238 + ^252 + ?7336, 

Vh = 2r749 + 2r7i47 + 277245 +^7343 • 

Here we have used the short notation 777V defined in Eq. with the omitted value M = 
M1M2M3 = 686. 

4.3 Numerical examples in 7-7-14 model 

Here, we give examples of numerical studies by using the 7-7-14 model, which is discussed in 
the previous subsection. For such studies, the numerical values of rjN defined in Eq. ([9]) are 
useful. The iV-dependence of 77 at is shown in Fig. [H 

We assume that both the up-sector and the down-sector of quarks as well as their Higgs 
fields have the Yukawa matrix, which is led in the 7-7-14 model. Such situation is realized in 
the case that we start with the U{8) gauge group and break it to ?7(4) x U{2)i x f/(2)^ by the 
magnetic fiux, and then the Pati-Salam gauge group is broken to the Standard gauge group 
by assuming VEVs of Higgs fields on fixed points. Alternatively, we break the f/(8) gauge 
group to f/(3) X U{l)i X U{2)l x U{1)2 x U{\)z by magnetic fiuxes and orbifold projections 
as discussed in section 3. Then, both the up-sector and down-sector of quarks can correspond 
to the Yukawa matrix led in the 7-7-14 model, although the up-sector and down-sector can 
generically correspond to different patterns of Yukawa matrices. In both cases, VEVs of the 
up-sector and down-sector Higgs fields are independent. 

First, we consider the case that VEVs of H^, and if^ are non- vanishing and the other 
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VEVs vanish. In this case, the relevant Yukawa couphngs are 

-Vc \ 



-Ve 



-T2y^Hl -yaHl 



\ 







Let us assume {H^) = —{Hi) for their VEVs. Then, quark mass ratios are obtained from these 
matrices as 



{my,,mc,mt)/mt 
{md,ms,mb)/mb 



(7.6 X 10-^ 6.8 X 10^1 1.0), 
(7.5 X 10-^ 5.1 X 10-2, 1.0), 



for r = i. Furthermore, the mixing angles are obtained as 

|Vci<:M| ~ 



0.97 0.24 0.0025 
0.24 0.95 0.20 
0.046 0.19 0.98 



Similarly, for r = 1.5i, quark mass ratios are obtained as 



{mu,mc,mt)/mt 
{mb,ms,md)/mb 

and the mixing angles are obtained as 



(2.1 X 10-^ 1.8 X 10-^ 1.0), 
(1.4 X 10'\ 1.7 X 10-^ 1.0), 



/ 0.99 0.13 0.00029 

0.13 0.98 0.13 
V 0.017 0.13 0.99 



Let us consider another type of VEVs. We assume that VEVs of H^, and are non- 
vanishing and the other VEVs vanish. Furthermore, we consider the case with {H^) — —{H^) 
and (H^) = (i?J)/3. In this case, the mass ratios are given by 

{mu,mc,mt)/mt ~ (2.9 x 10-^ 2.5 x 10"^ 1.0), 
{md,ms,mb)/mb ~ (4.4 x 10'^ 0.18, 1.0), 



for T = i, and the mixing angles are given by 



/ 0.98 0.22 0.018 
0.22 0.98 0.0014 
^ 0.017 0.0052 1.0 
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Similarly, for r = 1.5i the mass ratios and the mixing angles are given by 

{mu,mc,mt)/mt ~ (5.6 x 10"^ 4.7 x 10"^ 1.0), 
{md,ms,mb)/mb ~ (3.3 x 10"^ 7.1 x 10~^ 1.0), 



I VcKM I 



/ 0.98 0.22 0.0034 
0.22 0.98 0.000081 
V 0.0033 0.00081 1.0 



Thus, these values can realize experimental values of quark masses and mixing angles at a 
certain level by using a few parameters, i.e. r and a couple of VEVs of Higgs fields. If we 
consider more non-vanishing VEVs of Higgs fields, we could obtain more realistic values. For 
example, we assume that VEVs of H^, H^, H^, H} and Hj are non-vanishing and they satisfy 
-{H^) = {Hi) = {Hi) and {H}) = -{Hj)/2 while the other VEVs vanish. For r = 1.5z, we 
realize the mass ratios, niu/rnt ~ 2.7 x 10~^, nic/mt ~ 3.5 X 10-^ nid/mb ~ 7.3 x 10"^ and 
ms/rrib ~ 7.5 x 10~^, and mixing angles, \4s ~ 0.2, Vcb ~ 0.03 and Vub ~ 0.006. When we 
consider more non-vanishing VEVs of Higgs fields, it is possible to derive completely realistic 
values. Similarly, we can study other classes of models and they have a rich flavor structure. 



5 Conclusion 

We have studied three generation magnetized orbifold models. We have classified their flavor 
structures and studied explicitly possible patterns of Yukawa matrices. Our models have a 
rich flavor structure, especially compared with the corresponding models without orbifolding. 
Realistic quark masses and mixing angles can be derived within the framework of magnetized 
orbifold models. We can extend our numerical studies including the lepton sector. 

Here, we have studied the models, where all of three generations are originated from bulk 
modes. However, we have degree of freedom to put some of three generations of quarks and 
leptons on certain orbifold fixed points. In addition, we can assume that some Higgs fields are 
localized on certain orbifold fixed points. In such cases, we would have more variety of flavor 
structure. Furthermore, it is possible to consider localized magnetic fluxes on orbifold fixed 
points, which are independent of the bulk magnetic flux§ Since such localized magnetic fluxes 
would affect profiles of zero-modes, that is one of interesting extensions of our models. 

We have restricted ourselves to Abehan fluxes, but we can also extend our analysis to models 
with non-Abelian fluxes, which can reduce ranks of gauge groups. Moreover, although we have 
concentrated on the factorizable torus, (T^)^, it would be interesting to study possibilities for 
extensions to non-factorizable orbifolds [27] . 
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A Possible patterns of Yukawa matrices 

In this appendix, we show explicitly all of possible Yukawa matrices for 15 classes of models in 
Table [3] except the models with ij^l'^ = and the model without zero-modes for the Higgs fields. 

A.l (Even- Even- Even) wavefunctions 

Here, we study the patterns of Yukawa matrices in the models, where zero-modes of left, right- 
handed matter fields and Higgs fields correspond to even, even and even functions, respectively. 

A. 1.1 4-4-8 model 

Let us study the model with d/^b'*!) l-^m'*!; l-^bcl) — (^j 4, 8). The following table shows zero-mode 
wavefunctions of left, right-handed matter fields and Higgs fields. 














0U,4 


0U,4 


0U,8 


1 

2 


(01,4 + 03,4) 
02,4 


02,4 


(01,8 + 07,8) 

(e^'S + 96.8) 


3 






^ (03,8 + 05,8) 


4 






04,8 



This model has five zero-modes for the Higgs fields. Yukawa couplings Yij^LiRjEk are given by 

/ yaHo + y^Hi y^Hs + y^Hi y^H2 

YijkHk = y^Hs + ybHi -j^^y^ + ye)H2 + ydHo + H^) y^H-i + ydHi 
\ yA ybH^ + VdHi yeHo + yaH^ 

where 

ya = r]o + 27732 + r]64, yb = r]4 + r]28 + me + Veo, 

Vc = + + + Vd = Vl2 + V20 + V4.'i + V52, 

ye = 2?7i6 + 27748, 

in the short notation tjn defined in Eq. (Q with M = M1M2M3 = 128. 
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A. 1.2 4-5-9 model 



Here we show the model with |, \Icl^\, Ihc \) = (4,5,9). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 



ab\ 



00,4 
1_ (01,4 + 03,4) 
02,4 



V2 



00,5 
1 (01,5 + 04,5) 



V2 



^ (02,5 + 03,1 



00,9 

^ (01,9 + 08,9) 

^ (02,9 + 07,9) 

^ (93-9 + e6'9) 

1 (04,9 + 05,9) 



V2 



This model has five zero-modes for Higgs fields. Yukawa couplings Yij^LiRjUk are given by 



where 



/ ?7o v^r/ae \/27772 

^2^745 ^79 + ^781 V27 + ^?63 
V ?790 ^2^754 a/2?7i8 



/ 75(^720 + ^40) 



+ ^776 



^?32 + ^?68 



V 



^75 + ^785 

72^750 



75(??31 + ?741 + ^749 + ^759) 73(^13 + ^723 + ^767 + ^?77) 



^744 + ^764 
7744 + ^764 



V22 + V58 
VS + ^728 



V 

/ 779(^60 + ^80) 



4 



75(^71 + Vi9 + Vn + ^789) 71(^7 + ^737 + ^753 + ^73) 

^726 + ^746 V62 + ^82 

V24 + ^784 ^712 + ^748 

75(^721 + 1139 + r]51 + ^769) 73(^73 + ^733 + ^757 + ^?87) 

Ve + ^726 ^742 + ^778 

^716 + V5e V52 + ^?88 

73(^711 + V29 + V61 + ^779) 73(^77 + ^743 + ^747 + ^83) 

?734 + ?774 m + ^38 

in the short notation r]N defined in Eq. iQ with M = M1M2M3 = 180. 



V 



^715 + ^775 

^2^730 



^ 73(??60 + ^80) 

^725 + ^765 
V ^2^70 



A.1.3 4-5-1 model 

Here we show the model with (|/afe''|, {Ica'l, \Ibc\) = (4,5,1). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs field. 














00,4 


00,5 


00,1 


1 

2 


(01,4 + 03,4) 
02,4 


(01,5 + 04,5) 
^ (02,5 + 03,5) 
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This model has a single zero-modes for the Higgs field. Yukawa couplings YijkLiRjHk are given 

/ yo a/2?74 y/2r]8 \ 
YijkHk = V2r]5 + r]g) (773 + 777) Hq. 
\ ?7io V^Ve V2r]2 J 

Here we have used the short notation rj^ defined in Eq. (Q with the omitted value M = 
M1M2M3 = 20. 

A. 1.4 5-5-10 model 

Here we show the model with d/^b'*!, \Icl^\, l-^bPl) = (5,5, 10). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 














90,5 


90,5 


00,10 


1 




^ (61'^ + 


^ (01,10 + 09,10) 


2 


^ (02,5 + 03,5) 




(02,10 + 08,10) 


3 


^ (03,10 + 07,10) 


4 






_^ (04,10 + 06,10) 


5 






05,10 



This model has six zero-modes for Higgs fields. Yukawa couplings YijkLiRjHk are obtained as 

= ybHi + yeH^ ycHo + ^{yaH2 + y/Hs) + yaH^ -^{yaHi + yeH2 + ybH^. + ycHi) 
y ycH2 + ydH'i -^{ydHi + yeH2 + ybH^ + ycH^) ybHo + -^{yjHi + yaH^) + yaH^ j 

ya = r]0 + 2?750 + 2?7ioo, Vb = m+ ^745 + ^755 + ^795 + ^?105, 

yc = ?7io + ^?4o + ^?60 + '790 + '7iio, yd = r]i5 + V35 + V65 + Vsb + V115, 

Ve = ?720 + ^?30+^?70 + '780 + '7l20, Vf = '^V25 + '^Vlb + Vl25, 

in the short notation r]N defined in Eq. with M = M1M2M3 = 250. 
A. 2 (Even-Odd-Odd) wavefunctions 

Here, we study the patterns of Yukawa matrices in the models, where zero-modes of left, right- 
handed matter fields and Higgs fields correspond to even, odd and odd functions, respectively. 

A.2.1 4-7-11 model 

Here we show the model with \Ibc\) = (4, 7, 11). The following table shows zero- 

mode wavefunctions of left, right-handed matter fields and Higgs fields. 
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ab\ 



00,4 
02,4 



06,7) 
05,7) 
04,7) 



bc\ 



9 



10,ll^ 



^ (02,11 _ 09,11) 
(03,11 _ 08,11) 

^ (e^'ii - 9^'": 

_1= (05,11 _ 06,11) 



This model has five zero-modes for the Higgs fields. Yukawa couplings Yij^LiRjUk are given by 



where 

72(^74 - ?7l36) 
^?81 - - V95 + V73 
y2('7i50 - ^?18) 

V2(^?80 - r]52) 
VS - r]25 - 11129 + ^7151 
\/2(^?74 - ^7102) 

/ \/2(?7i44 - ^732) 
^787 - ^109 - Vi5 + V67 

^2(^710 - ^122) 

^2(^148 - ^7104) 
^7171 - ^115 - ^39 + ^17 
V2(^94 - ^?38) 

^2(^724 - ?7l08) 
^753 - ^?31 - ^123 + ^?101 

V2(^?i30 - me) 
in the short notation tin defined in Eq. (| 



ylj 



V2{7]92 
V139 - V29 - 

V2(r?62 

\/2(?7l46 
72(^76 

m - viu - 

V^(^78 

V2(^?148 
^783 - ^727 - 

Vl41 - V57 

72(^790 
91) with M -- 



- ^48) 
^7125 + Vlb 

-Vie) 

- rjii + ^769 

- ^?118) 

- ^7120) 

^743 + ^153 
-^34) 

- ^7104) 
' ^7127 + ^771 

- ^50) 

- ^20) 

- V134) 
= M1M2M3 



V2{Vl28 


-^40) 


V51 - Vll7 - 


- ^737 + ^7103 


V2iV26 


- ^114) 




- ^7124) 


Vl35 - Vl07 


- ^?47 + ^?19 


72(^58 


-^30) 


V2{vi2 - 


- Vloo) 


V89 - V23 - 


^7131 + V65 


\/2('7l42 




V2{V72 


-Via) 




Vl3 + Vli9 




- Visa) 


V^{Vl52 


- V68) \ 


V79 - ^145 - 


- V9 + V75 


V2iv2- 


- V86) J 


= 308. 





A.2.2 4-7-3 model 

Here we show the model with (l/afe"*!) l-^ca"*!? l-^bcl) — (4,7,3). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 





L,i\-') 




Hk{\'') 




1 

2 


00,4 
02,4 


_^(0l,V_06,V) 

^ (93-5 - 9^'^) 


^(01,3 _ 02,3) 



This model has a single zero-modes for Higgs fields. Yukawa couplings YijkLiRjHk are obtained 
as 

2 / 72(^74 - ^32) V2(^?2o - Vs) a/2(^?4o - Vw) \ 

Yi'^Hk = -Tk^oI V17 + V25 - Vii - V31 Vl + Vii - Vi3 - V29 Vi9 + V23 -V5-V37 , 

^ \ V2{1]38 - Vio) V2{V22 - V3l) \^{V2 - V26) / 
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in the short notation 7]^ defined in Eq. with M = M1M2M3 = 84. 



A. 2.3 4-8-12 model 

Here we show the model with \Icl^\, 1^11^1) = (4,8, 12). The following table shows zero- 

mode wavefunctions of left, right-handed matter fields and Higgs fields. 



00,4 
92,4 



V2 



6,8\ 



73 - 

_1= (03,8 _ 05,8 



_1= (92,12 _ 010,12 
_1= (03.12 _ 09,12 
_1= (04.12 _ 08,12 
(05,12 _ 07,12 



This model has five zero-modes for the Higgs fields. Yukawa couplings YijkLiRjHk are given by 



where 



71 



















yd 


-1// 


Vb 



71(2/6 - y/i) 







73(^/-^') 

V 



Vc - yk 







Vc - Vk 









-^{vb-yh) 



\ 



J 



and 



Va 


= ^0 +^96 +^7192 +^?96, 


Vb 


= ^74 + ^7100 + ^188 +^792, 


Vc 


= ^8+^04+^7184+^88, 


Vd 


= ^712 + ^7108 + ^180 + ^784, 


Ve 


= ?7i6 + ^112 + rjn& + ^780, 


Vf 


= ^?20 + ^?116 + ^7172 + ^?76, 


Va 


= ?724 + ^?120 + ^168 + ^772, 


Vh 


= ^?28 + ^?124 + ?7l64 + VG8, 


Vi 


= ?732 + ??128 + ^?160 + ^764, 


Vj 


= V36 + ^7132 + ?7l56 + Vm, 


Vk 


= '740 + ^136+^?152 + '756, 


yi 


= V44 + VmO + VUS + V52, 


Vm 


= ?748 + 77144 r7l44 + ^748, 







m 



the short notation r]N defined in Eq. with M = M1M2M3 = 384. 



A.2.4 4-8-4 model 



Here we show the model with (1/^6 |, \Icl^\, Ihc \) — (4,8,4). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 
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ab\ 



00,4 
02,4 



6,7\ 



^ (03,7 _ 05.7 



5(0 



bc\ 



1,4 



e 



3,4^ 



This model has a single zero-modes for Higgs fields. Yukawa couplings YijkLiRjHk are obtained 
as 



Hn 



[ Vb -Vc 

^^{ya-vd) 

V -2/c 



where 



Vc 



VO + 2l]32 + VGA, Vb 
VU + V20 + V'l^ + V52, Vd 



^?4 + ?728 + '736 + ??60, 

2r]iQ + 27748, 



in the short notation r]N defined in Eq. with M = M1M2M3 = 128. 



A.2.5 5-7-12 model 

Here we show the model with (| |, {Icl^, 



^bc\) = (5, 7, 12). The following table shows zero- 



mode wavefunctions of left, right-handed matter fields and Higgs fields. 







RjiX"") 


HkiX'^) 





00,5 


M^'' - 


06,7) 


^(01,12 _ 011,12) 


1 






05,7) 


_1. (02,12 _ 010,12) 


2 


^ (02,5 + 03,5) 




04,7) 


(03,12 _ 09,12) 


3 




_1= (04,12 _ 08,12) 


4 








_1= (05,12 _ 07,12) 



This model has five zero-modes for the Higgs fields. Yukawa coupling YijkLiRjHk are given by 



^ ijk-'^-^ k 



where 



1 
71 



/ A/2(r75 - 7765) ^2(^185 - ^7115) \/2('755 + ^7125) 

^173 - ^103 - ^187 + ^163 ^67 " ^?137 " ^53 + ^17 ^7113 - ^43 " ??127 + ^197 
\ ^779 - ^?149 - ^19 + ^789 ^101 " %1 " ^199 + ^151 ^7139 - ^7209 " VaI + ^?29 

I ( V^iVno - Vno) v^(?7io - Viso) \/^{Vi90 + V5o) 

V2- Vl42 - V58 + V82 ^178 " ^738 " ^7122 + ^7158 ^762 " ^?202 " VU8 + V22 
\ Vwe - V26 - Vl94 + V94: V74 " ^7206 " ^746 + ^794 ^7106 " ^734 " ^7134 + ^7146 



1 

71 



^ 72(^775 - ^?135) 72(^7165 - ^?45) V2(?7l5 - ^7195) 

^177 - ^?33 - ^7117 + ^?93 V3 " ^?207 " ^?123 + ^787 ^183 " ^27 " ^?57 + ^7153 

\ V9- V201 - V51 + Vsi Vni - V39 - Vi29 + Vsi Vm - Vui - Vui + V99 
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V2('7ioo - Vmo) 

-7^ \ V68- V2O8 - ^128 + ^7152 Vn2 

^2(^7145 - to) 
—7^ I ^107 - ^37 - ^747 + ^23 ^73 - ^143 - ^193 + ^7157 ^167 



V^iVso - V200) \/2(r7i6o - V20) 

V32 - V52 + V88 V8 - ^?148 " ^7188 + ^792 

^7116 + ^164 Vne - Vi04 - + Vre 

V2 (r^gs -V25) (^85 - ^?155 ) 

^97 - Vl3 + ^783 
?7l81 + ^169 



in the short notation t]n defined in Eq. with M = M1M2M3 = 420. 



A. 2.6 5-8-13 model 



Here we show the model with (|/^^''' 



'ab h l^^a^l. \Hc \) = (5,8, 13). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 



ab\ 



90,5 



V2 



06,8) 
95,8) 



bc\ 



f 

V2 
1 

V2 



V2 
1 

V2 



(01.13 . 


_ 012,13 


(02,13 


_ 011,13 


(03,13 


_ 010,13 


(04,13 


- 99.13^ 


(05,13 


_ 08,13^ 


(06,13 





This model has six zero-modes for the Higgs fields. Yukawa couplings YijkLiRjHj. are given by 
Yt^H, = ylH, + y]^H,+ylH, + ylH, + yt,H, + ylH,, 

where 



1 
1 

71 



^ \/2(?75 - 77125) 

^?203 + ^213 - ??83 - ^187 
\ ^109 - ^21 + ^799 - ^7229 

/ v^lto - ^775) 

% + ^7237 - ^7133 - ^7107 
\ ^?211 - ^179 + ^101 - ^?29 

/ y2('7ii5 - to) 

^7197 + ^137 - ^67 - ^793 
V ^11 - ^7141 + ^219 - ^171 

V^('785 - ^745) 
^123 + ^163 - ^253 " ??227 
^189 - ^?59 + ^19 - ^?149 

V2{r]235 - ^155) 

^77 + ^157 - ^53 - ^27 
^7131 - ^7259 + ^181 - V51 



v^l^/wo - ^770) 

^7122 - ^7138 + ^?18 - ^242 
^786 - ^174 + ^7226 " ^734 

v^('7io - to) 

^198 - ^62 + ^132 - ^152 
?7ll4 - ?7l46 + ^94 - ^?166 

^2(^7210 - ^750) 
^72 - ^?258 + ^7102 - ^7158 
^206 - ^754 + ^?106 - ^7154 

V^l^/llO - ^7150) 
^202 - V58 + ^?98 - ^162 
^76 - ^254 + ^214 - ^746 

V^('790 - ^7170) 
^7118 - ^7142 + ^?222 - ^?38 
^194 - Vm + Vl4 - V246 



V73 - V57 + Vm - V47 
V239 - Vl5l + V31 - VWI 

V2{Vl85 - V55) 
V127 - V257 + V23 - V153 
VSI - V49 + V23I - ^7159 

V^ivib - V145) 

V193 - V63 + V223 - Vl67 

Vu9 - V249 + V89 - mi 

V2{vi35 - to) 

^77 - ??137 + ^797 - ^733 
^201 - ^771 + ^111 - ^241 

V2(^?i05 - ^25) 

^7207 - ^7183 + ^?103 - ^?233 
^71 - ^129 + ^7209 - V79 
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I ( V2iV35-Vi65) \/2(ry23o - ^730) V2(^?95-te) \ 

^ 1^69 - Vei + r/139 - V2,i Vi2e - Vi3. + Vise - Vr. - + V9 - V121 J 
in the short notation r]N defined in Eq. with M = M1M2M3 = 520. 

A.2.7 5-8-3 model 

Here we show the model with \Icl^\, l^bcl) = (5,8,3). The following table shows zero- 

mode wavefunctions of left, right-handed matter fields and Higgs fields. 













1 

2 


90,5 

^ + e''') 


^(e3._ 


97,8) 
96,8) 
95,8) 


^(91,3 _ 9^,3) 



This model has a single zero-mode for the Higgs field. Yukawa couplings Yij^LiRjEj. are given 
by 

I ( V^ijlh - ^35) V2('750 - ^710) y2('725 - ^755) \ 

^ij^k = -J=\ VA3- V37 - Vl3 + ^753 ^72 " ^738 " ^58 + ^722 ^747 -^77- Vu + V23 , 
\ ^729 - Vll - V59 + Vl9 VA6 - V3A " Vm + V26 Vl " ^741 " ^31 + ^749 / 

in the short notation r]N defined in Eq. with M = M1M2M3 = 120. 
A. 3 (Odd-Odd-Even) wavefunctions 

Here, we study the patterns of Yukawa matrices in the models, where zero-modes of left, right- 
handed matter fields and Higgs fields correspond to odd, odd and even functions, respectively. 

A. 3.1 7-7-14 model 

Here we show the mo del with = (7,7,14). This model is studied in the 

subsections 14.21 and 14.31 in detail. The zero-mode wavefunctions of left, right-handed matter 
fields and Higgs fields are shown in Table HI 

This model has eight zero-modes for the Higgs fields. Yukawa couplings YijkLjRjHk are 
obtained as 

Yi.kHk = vIHq + yl^Hi + ylH2 + y^H^ + yt^H^ + y^H^ + yJiYe + 4^^7, 
where yf- is shown in Eq. <^ with M = M1M2M3 = 686. 
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A. 3.2 7-8-15 model 



Here we show the model with (l/ab"*], \Icl^\, I^Pl) = (7, 8, 15). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 



(02,7 _ 05,7) 



(01,8 _ e7,«) 

_^ (02,8 _ 06,8) 
1 (03,7 _ 05,8) 



V2 



00, 



15 



_^ (01,15 + 014,15) 
_^ (02,15 + 013,15) 
(03,15 + 012,15) 
^ (04,15 + 011,15) 
^ (05,15 + 010,15) 
^ (06,15 + 09,15) 
^ (07,15 + 08,15) 



This model has eight zero-modes for the Higgs fields. Yukawa couphngs YijkLiRjHk are given 

by 



where 




' r\22h - r)l5 ^330 - ^90 V405 " ^7195 
^7345 - ^7135 ^7390 " ^730 ^285 " V75 
\ V375 - V255 V270 " Vl50 VW5 " »745 

??113 - V97 - Vl27 + V337 V218 " ^7202 " ^722 + %98 %23 " %07 " 

?7233 - ''723 - ?7247 + ?7383 ?7338 " ?782 " ?7l42 + ?7278 ?7397 " ?7l87 " 

?7353 - ?7l43 - ?7367 + ?7263 ?7382 " ?738 " ?7262 + ?7l58 ?7277 " ?767 - 

f Vl- V209 - r)239 + ^391 ^106 " ^314 " ?7l34 + r)286 ^211 " ?7419 " 

?7l21 - ?789 - ^7359 + ^7271 ^7226 " ^7194 - ^7254 + ?7l66 ^7331 - ^7299 " 

\ rj241 - r}31 - ^361 + ^151 ?7346 - ?774 - ?7374 + V^e ?7389 - ?7l79 " 

/ ^7111 - 77321 - 77351 + 7/279 Vd " ^414 - ^7246 + ^174 ^799 - 7/399 " 

V9 - V2OI - ^7369 + ^159 ^114 - ^7306 - ^7366 + ^754 ?7219 " ^411 

\ ^129 - ^81 - ^249 + ^39 ^234 " ^186 " ^354 + %6 ^339 " ^291 " 

?7223 - ?7407 - V377 + Vl67 VU8 " %02 " %58 + V62 Vl3 " ^197 

?7l03 - ?7313 - ?7257 + ?747 ?72 " ?7418 " ?7362 + ?758 ?7l07 " ?7317 

?7l7 - ?7l93 - ?7l37 + ?773 ?7l22 " ?7298 " ?7242 + ?7l78 ?7227 " ?7403 

/ ?7335 - ?7295 - ?7265 + ?755 ?7230 " ?7l90 " ?7370 + ?750 ?7l25 " ?785 - 

^215 - ^415 - ^145 + ^65 VUO - V310 " ^7250 + ^7170 ^5 " ^205 - 

V ^95 - ^305 - ^725 + ^185 ^10 " ^410 ^ ^130 + ^290 ^115 " ^325 " 

^393 - ^183 - ^153 + ^57 ^342 " ^778 " ^258 + ^162 ^237 - ??27 " 

^7327 - ?7303 - ^733 + ?7l77 ?7222 - ^7198 - ^138 + ^7282 ^7117 - 7793 - 

?7207 - ?7417 - ?787 + ?7297 ?7l02 " ?7318 " ?7l8 + ?7402 ?73 " ?7213 " 



1 

1 
72 



- V83 + V293 

- ?737 + ?7l73 

- ?7l57 + ?753 

-^729 + ^181 ^ 

- ?7l49 + ^761 

- ?7269 + ^59 / 

- ^141 + ^69 

- ^7261 + ^751 

- ^381 + ^71 

- ^7253 + ^43 \ 

- ?7373 + ?7l63 

- ?7347 + ?7283 / 

- ?7365 + ?7l55 ^ 

■ ^355 + ?7275 

- ^235 + ^395 / 

?7363 + ^267 
^7243 + 7/387 
?7l23 + ?7333 
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I ( V281 - Vn - V41 + Vim Vsse - %4 - Vug + V274 V349 - V139 - V251 + ^7379 \ 

= ~in ^401 - Vvn - V79 + V289 ^?334 - VS6 " ^26 + ^7394 V229 " Vl9 " ^7131 + VSH 

\ V319 - V3U - Vl99 + Vm9 V21i - V206 ~ V94 + ^326 ^109 - VWI - Vll + V221 J 

in the short notation r]N defined in Eq. with M = M1M2M3 = 840. 



A. 3.3 7-8-1 model 



Here we show the model with (1/^6 |, {Icl^, Ihc \) — C'')8,l). The following table shows zero- 
mode wavefunctions of left, right-handed matter fields and Higgs fields. 



ab\ 



-^(0178 _ 07,8) 
(03-8 _ 05,8) 



bc\ 







0,1 



This model has a single zero-mode for the Higgs field. Yukawa couplings YijkLiRjHk are given 

by 



' V2{r]5 - 7735) 

^743 - ^737 - ^713 + ^753 V2 " 
V ?729 - ^711 - ^759 + ^719 



V2iV50 - Vw) 



V2('725 - ^755) 



^738 - ^58 + ^722 ^747 -^77" Vl7 + ^723 
- ^34 - Vl4. + V26 Vl - VAI - ^31 + ^749 



in the short notation r^^r defined in Eq. ([9]) with M = M1M2M3 = 56. 



A. 3.4 8-8-16 model 



Here we show the model with \Ica^\, \Ibc\) = (8,8, 16). The following table shows zero- 

mode wavefunctions of left, right-handed matter fields and Higgs fields. 




1 
2 
3 
4 
5 
6 
7 
8 



^ (02,8 _ 06,8) 
_1= (03,8 _ 05,8) 



WW 



V2 



t 

x/2 



J= (03.8 



07,8) 
06,8) 
05,8) 



bc\ 



00,16 
(01,16 + 015,16) 
J= (02,16 + 014,16) 



V2 

72 + 

72 + 
72 + 

_1^ (06,16 + 
^(07,16 + 
08,16 



013,16) 
012,16) 
011,16) 
010,16) 
09,16) 



This model has eight zero-modes for the Higgs fields. Yukawa couplings Yijj^LiRjHk are obtained 



as 



Y,';Hk = ylH, + ylH, + ylH2 + ylH, + yt,H, + yl^H, + ylH, + ylHT + ylH, 
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where 



and 



4 







-Vf 





-Vf 





-yd 





-yd 






ya 
yt 



yc 
yd 
ye 
yf 

yg = 
yh = 

yi = 



= r]o + 2(r7i28 + 2?7256 + 2?7384) + V512, 

= V8 + V12O + Vl36 + ^248 + ^7264 + ^376 + ^7392 + ^504, 

= ^716 +11112 + VU4 + V240 + r]272 + ?7368 + V^OO + ^7496, 

= ^724 + VWi + " 



^7156 + ^?232 + ^7280 + ^7360 + ^408 + ^7488, 
V32 + V96 + VW4 + V224 + V288 + V352 + V416 + ^480, 



= ^740 + ^?88 + + r]216 + r]296 + V344 + V424 + ^472, 

^748 + ^?80 + 



^180 + ^7208 + ^?304 + ^7336 + ^?432 + ^7464, 
V56 + V72 + VlSS + V200, +V312 + V328 + 



2(??64 + ?7l92 + ?7320 + ^?448), 

in the short notation rj^ defined in Eq. iQ with M = M] 



?7440 + ^456, 

M2M3 = 1024. 
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